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SECTION I 


INTRODUCTION 


A method of analytically determining various foil bearing 


characteristics for the foil thrust bearing and the foil journal 


bearing were developed during this program. The approach used 


in developing the analytical programs is discussed in Volume I. 


This volume presents the three computer programs in Fortran 


source listing form and a User's Manual for each. 


The 


III, and 


(a) 


(b) 


(c) 


= oat 
¥ ° 


three programs, discussed separately in Sections II, 


IV are: 

Foil Journal Bearing Elasticity Computer Program 
"Foils". 

Foil Journal Bearing Elasticity/Hydrodynamic Computer 
Program 

Foil Thrust Bearing Elasticity/Hydrodynamics Computer 


Program 
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SECTION II 


FOIL JOURNAL BEARING ELASTICITY COMPUTER PROGRAM "FOILS" 


2.1 General 


The efforts to analyze foil journal bearings required the 
formation of a specialized nonlinear elasticity solution. Program 
"Foils" has been developed and forms the basis of a coupled 
elasto~hydrodynamics solution to foil journal bearings. Rotor 
dynamics simulation requires accurate prediction of bearing stiff- 
ness characteristics; therefore, this nonlinear elasticity solu- 


tion is presented as a separate program. 


Program "Foils" is a static elasticity solution for foil 
journal bearings. No inertia or hydrodynamics effects are incor- 


porated. However, testing of the foil bearing test rig has demon- 


strated tha! location of critical speeds can be accurate. - 
predicted indamped rotor dynamics model utilizing bearing 
stiffnes ed by Program “Foils”. 


2.2 Program Limitations and Assumptions 


6, 8 or 10 nodes per foil 

A maximum of 160 nodes 

The coating is assumed to be nonstructural 
Constant thickness 

Misalignment not permitted 


No frictional effects considered 


ooo GO 0 OO OC 


The shaft and bearing housing are assumed to be rigid 
2.3 Theoretical Basis and Flow Chart 
Some typical foil bearing geometries are illustrated in 


Figure 1. The foil elasticity is represented as a series of 


short curved beam elements with the coordinate system and sign 


convention of Figure 2. 
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2. Foil Geometry Coordinate System. 
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: 


By application of Castigliano's second theorem, the element 
and system stiffness matrices may be determined. The solution to 
the stiffness equation is provided by a straight-forward displace- 
ment matrix approach. The attachment of the foil to the bearing 
housing is assumed to be a “pinned joint". The difficulty in 
the solution is encountered when the foils are loaded. The 
foils begin to contact each other and encounter kinematic con- 
straints of the journal and bearing housing. It is apparent that 
this provides a highly nonlinear stiffness characteristic of foil 


journal bearings. 


A simplified flow chart of the geometry/elasticity program 
"Foils" is presented in Figure 3. The following is a discussion 
of this program and should be useful in understanding how the pro- 
gram deals with foil overlap and boundary conditions. A Fortran 
listing is provided in Appendix A. 


2.4 Assembly "Geometry Routine" 


As all foil displacements must be referenced from a pre- 
scribed datum, the assembly routine provides such a datum; the 
unloaded assembly position. The predicted assembly position rep- 
resents the geometry of the foil bearing prior to shaft insertion. 
This geometry routine has been programmed to search for an accept- 
able, constant curvature position, Figure 4. When Ro (radius of 
curvature at the assembly position) satisfies the constraint that 
R' = R+t where t is the foil thickness, the assembly subroutine 
reports Ry as well as Ry (the radius to each node from the geo- 
metric center of the bearing). It should be noted that Ro may or 
may not be equal to the free foil radius. (Generally the larger 
the number of foils, the more likely Ry will be the free foil 


radius.) 


2.5 Elasticity "Routine" 


As noted in the flow chart (Figure 3), the elasticity solu- 


tion commences by obtaining the symmetric "preload" position of 
5 
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NODE N 
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the foils. This is the configuration when the shaft is inserted 
in a zero-"g" position. The program proceeds from there to obtain 


the eccentric position of the foils while checking for: 


(a) Constraint or release of nodes on adjacent foils 


and/or contact with boundaries 
(bo) Constraint or release of nodes on the shaft (S) 


(c) Constraint or release of nodes on the bearing 


housing (B) 


Nodal constraint is provided when displacements indicate a viola- 
tion of kinematic boundaries. Release of nodes occurs when the 


traction loads produced by constraint are not physically realizable. 
2.6 Data Presentation 


The geometry and coordinate system are shown in Figure 5. The 
node number system for a typical 8 foil geometry with 10 nodes per 
foil is illustrated in Figure 6. The numbering system is automatic 
Once the number of nodes per foil and number of foils have been 


specified. 


Figures 7 through 13 exemplify the output or data presentation 
of Program "Foils". Figure 7 is a summary of the input data as 
well as the specification of the assembly position, radius of cur- 
vature (Ros Figure 4) and the nodal radii (Rye Figure 4). In addi- 
tion, the shaft radius which will be utilized to determine the sym- 
metric preload position as well as the eccentricity control are 
presented. "PSI" is the angle (¥) of Figure 5. Figures 8 through 
10 present the results for ecc = 0.0 whereas Figures 11 through 13 
present the results for ecc = 0.002 in. The summaries of Figures 
10 and 13 note the nodal radial (total) load acting on the shaft 
due to the preload and eccentricity. In addition, the X and Y 
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Figure 5. Geometry and Coordinate System for Elasticity 
Program 
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Figure 6. Node Numbering 


10 62500 


Azeuung ezeq ynduy ‘4 eanbty 


S 
ae 
tS 
2 
“0 20-300002° 20-3000U2" 1o-3000EI° = TSd*LIN13°4d73U" IWNIS3 nN 
00se*2 = siIGvd L4vHs i 
errr rrr rrrrrrrrrrrrrrrrrTrirtiirrtirtiiricooco rl ke eee eee ee 
UT 
6 
4 
] 
9 
S 
f - 
£ 
2 a 
Owi2°2e t 
rer rrrrrrrererrrerrrrrrrrrir irri titi tilili ir he. kk ake 
ONT Hya8 40 H34N39D WOHS SNLOVH A laWSsSS¥ 4QON 


T0*300t%° = NOILISOd ATHARSSY ap SHfilyAuIID 40 SNTQVYa 


*NOTLISOd ATAWSSSY 3HL HOS SANTA ONTMONIOS AHL SFAIN AGLaWUAD IIOS 
990 ANONT NYA *£3SVD GAVSBANUD vw SL SIHL 


AMLAWORD ATAW3SSY YNIGvad TOS 


peennceneneccsenoe cease eeeee SON neOSOEHHONHEH OOO ED PLODRERDO OLESEN HDOODAHDESHNHEDIEHHDHODINSOMANHHORAANHLOAHAHOBGELHLOSSSRLSOS 


zor azese" =1104 80d ‘3 
BU 3I0000E° = ALIDILSV 1a 30 SHOW 


PPPPPPTVTTETTPTTTerCrrrrrrrrrrrriorer err rie de ele 


SHHHOHOHOOHOSOOHEOHEOS 


000% *y 000%2°2 onTOU* oo2to* 00000°9 00167"? ul el 
4 VIOs 3384 HLONZIDNY T104 SSINMDIHL ONIL VOD SS3INWOIHL 104 HLONST TwlXy SNTQVM ONT MY as VWOs7S3NVUN 21103 30 °UN 


HOCHHOSOS HOV OOOOD 


POO REO OAH O OND OOO DU TIT DNDN HOON NOOO HOO OFODEOADO OE OOFOODUDEUODHIOGHHHOSHLHIHIFHVOOLOOEDS 
asai i104 21 


HONHHDOHEAHOES eonoo 


SYEGESTESSESTESGSESTSSTSESECSSSSFESLSEGTHSSHEITHGED ENS SESS SHSSFISSSTF HOGG S ASST SIIFIGS 


A ETTORE LL IN 


0°O = Aqtotazquecog AOF SsATNseyu °8 eanbtg 


vu7oe el” yae « uy te Pe | suue 
00v0" RIGe"e yl62°- pose) 3UU4 
0o0g"* oL52°2 ttte*- 64eu"2 3U0N 
vSll"2- 21y2°2 L601°= rt Gated 4u0N 
2%es" ot OvLer2 0000*%u Oric’Z JuON 
Sulvg dv 1H3A0 30uud IWwI0ve VN T sa 14390 VWwlOvn WILINIe = Vos 


VSL "EL aere*2 ROE *=- 2cR8"T 3U0N 
udiy "A= o0r2°2 81 OE°- 2hre y] 3UUN 
€ov0'= R2y2°2 £60%°= Steel 3UON 
oSlt"2 7B%2°2 980%" 96EB"T 4UON 
6l00°> 6€S2°2 198E°- LLGB*T SUUN 
eSeL°eI7 99s2°2 RQIE*= 2566°T 3U0N 
ovv0"> R9S2°2 ol62°- 8656" 3UON 
0000° nls2°2 {112°= 6S%u"? SUON 
69°27 2192°2 L60i1°> Sist’2 3UON 
2RTS*OT OL2°2 0000°O OvL2°2 4UON 
SHIvd dy T43AV 3odus Twldve NITY 1430 Widva AWTLINIG = lugs 


- 
wuwnag 
mene 


eomoowoaoonw 
- 


-anre 
mere 


7 


7SBL°EL Nev2°2 BU9E*= Z2ces't 
6994°8= 00v2°2 BI6E°= ZHoweT 
1000° Reve°2 €60%°~ Sees" 
ved"? 2By2*2 9R0%"= QnEee"t 
800° 6ES2°2 198 °= Liget 
L282" E17 09S2°2 g09e"= 2c6e'T 
00v0"= R952°2 ol62°- ReS6"] 
0000° alge°2 1G Ge didi 6S60°? 
969T°2- zi92°2 Lo01"- Sist*2 4aUON 

BBIS°OT OvLe°2 0000°v Nv.et2 3UON 

SulVd dv IH3AU 39du4d WIV" TnI dy 1430 Widva WILINIY = 1103 


3U0N 
dU0N 
saUON 
ZUU0N 
aU0N 
3U0N 
3UON 
SU0ON 


- 


omoomoeroot 
- 
nn 


" 


-“Qeauw €eagoe 
RAANNAAKRAAAH 
oon 


OLBLEL ofy2°2 BuIE*= 2zeRt 3auu’% 
vely R= C0n2°2 816E°- 2hos 3408 
T000" P2n2°2 €604°= Stes 3U0N 
9691 °2 28e2°e 980%"= QbEBTT au0% 
6100°- 6€S2°2 198E°= Lise ty 3UUN 
ursl "etm 0962°2 809E "= 2s6e'T 3UUN 
0000" RgG2°2 ol62°- ROSb°T 400 
0000" ols2°2 tlt2° 659U "2 3UON 
(99 2- z192°2 L601°- sist’e2 3UON 
eRtsat NeL2°2 000070 O+er? 43U0N 
sul¥g dV usAu 3ddu4 Wide WNT 4y 1430 TWw1dveg AWILINIe = Mod 


- 


uivd 

TWGON LOWLNOO 
INAGNadad UNV 
LNACGNaddONT 


- 


£ 
0 
0 
S 
0 
2 
i) 
0 
€ 
0 


uPBlL et Nene’? BO9ET = Zere I 
8999 °R= n0%2°2 BL6E°= anno" 
vo00" Revnz*2 £60%"= Steer 
7912 2Be2°2 Qa’ = gFee ) 
0000 '7 feo2'2 198E°= L198°] 
29RL"F17 09s2°2 809F*= Zeek 
0ou0” R992°2 0462°- ROS6"] 
0000" olg2°2 rie = 6Seu"2 
ONISNOH 6E9L 27 7192°2 L6Q,°> Sict’2 A3UON 
ONIUVad — 4 E9TS"Ot Oo1z°2 o000°v Oreer2 4UON 
HLIM *9°a © SulVq gv TH3aA0 3duus WAVY Wyt da 1430 WIOve WILINIn = VWg9 


LOWLNOO vooute = AL19101NI993 


3UON 
3UON 
3UUN 
SUON 
3u0N 
3UUN 
SUON 
3U0N 


LavHs 
HLIM 
LOWLNOO 


_ 


owe e ane 


" 


omoonomoon 
Oa 


-Rer e@auvarnagce 


NOLIalaysad W4ALSAS 


d$4i 104 21 


8 
*9*@ 


yn 


°9°a 


*9°9 


*9°8 


0°O = AQtoTtazqUSeDOq AOF sz[Nssy 


6 
0 


2 


SYIVd dV 1N3AU 


coro acoso 


~ 
Nu 


0 


- 


“ 


SUIVd dv 1u3AU 


Sulvg dv lua 


omooro noo ®f SPomoomorooan 


nw 


n~ 


0 


iy 


n 


Sul¥g d¥Iu3A0 


Sulvd dvi 


omoc oO >Porcooaomoor 


nw 


Nu 


fe) 


~ 


9818" EL ntre’?2 
S61G°R- 00%2°2 
1000° Ren2°2 
z2s8t"2 28%¢°2 
z020° RESe*2 
z2o28 Fn 09s2°2 
0000" A9S2°2 
0000" oLs2°2 
vLel°2- z192°2 
BSTS°OT Ov,2°2 
304043 Twluvs JINT Au 
lez3"Ft OEve*2 
S925°R- 00%2°2 
1000° R242°2 
OLel°2 7B%2°e 
2020" PES2%e 
llz8"et- 09S2°2 
tovo"= R9IG2°2 
0000° o4s2°2 
OS6[°2- ztgz°2 
602S°01 OvL2°2 
304u4 WWI1dvE JWNIdd 
zlz8°el 0€v2°2 
2925°R= 00%72°2 
1000° Rav2°2 
uS6T "2 2892°2 
ssuo” REG2°2 
TluB"Et> n9S2°2 
0000° R9G2°2 
0000°- 0lS2°2 
ISL1"2- 2192°2 
0025°01 OvL2°2 
39404 WWI10Vy Wwntsu 
6O6L "ET 0f92°2 
B68y°R- 00%2°2 
£000°> R2n2°2 
6olte2 28%2°2 
SSv0° RES2°2 
196L°ET- 09S2°2 
0000" R952°2 
0000°- olg2°2 
Olt 2- 2192°2 
SLIS*Ot OvL2°2 
30403 WIdovy WNT dy 
296L"ET oey2"2 
2984" B- oove"2 
£000°> R2v2°2 
6lLt"2 28e2°2 
Ssu0" RES2*2 
BS6l*FI- 09s2°2 
10vg"- R9G2°2 
0000°° nLo2°2 
badt*a- 7192°2 
“BIGOT OvLete2 
3dd04d Tvldvy Twn Tt da 
ES6L°ET nev2°2 
vEBe B= o0%2°2 
€000°° Rave2°2 
2l.°2 282°? 
SS00° RESe*e2 


°6 amnhtaq 


BO9Ee = 
RL 6E°= 
£60%°= 
QAny’ = 
198E C= 
RO9E* = 
ol62°- 
(itz 
L601°- 
0000°u 


1430 Wvldve 


B09E°- 
BL6E°= 
€60%°= 
980%°- 
198E°= 
RO9E "= 
oL62°- 
1tt2°- 
260T°- 
0000°v 


1430 Wwldve 


Bu9IE*= 
816E° 
e609" 
980%°- 
198E°= 
809E° 
0L62°- 
ttte-= 
L601°" 
0000°0 


W430 Wwidve 


809E*= 
B16E°- 
€60%°= 
980% 
198E "= 
B09E’= 
0L62°- 
Ttte>- 
L6ot'- 
0000°0 


1430 WvI0Ve 


BO9E*= 
8 16E°= 
€60%°= 
9ROe 
198ET= 
B09E"= 
0L62°° 
rite: 
L60T°> 
0000°0 


1490 WI0ve 


BYu9E’= 
RL 6E°- 
E60%°- 
ORDe*= 
LIE = 


ZchH Ty 
Zhen] 
Seavey 
Gee rT 
lige] 
2566") 
ROCb*] 
6S9u°? 
Sisl’2 
0%l2°2 
AILINIS 


2ceery 
agus *t 
SFEB*] 
96ee°T 
Lige*y 
Z2c6e tT 
RESE*] 
65%U°2 
Sisl*2 
O%l2°2 
WWILINIY 


ZebR*y 
Zeer] 
Sees y 
goed Ty 
LL96"1 
2o8b°T 
B6Cb"T 
6S60°2 
sist’2 
ooler2 
WILINIn 


2cee"| 
ZkoBeT 
SEER °T 
Soeur] 
LLge° 1 
2556" 
BES6*y 
65%0°? 
sist*2 
OvL1e°2 
WILINI 


2Zesey 
erokey 
Seen*y 
QAENTT 
LLgb 1 
etSRet 
BESS 1 
Buru'z2 
Sist*? 
Overe2 
Wiltnle 


ud 


ennne 


c 


3UUMN 
4du0% 
3UU0% 
SUON 
3UON 
AUON 
3U0N 
3U0N 
SUON 
3UON 
= W043 


S3U0N 
3u0Ow 
3U0N 
300N 
4aU0N 
4UON 
3U0N 
3UON 
3UON 
S3UUN 
= Vied 


3U0N 
3Uu0N 
3U0N 
3u0'y 
3UUN 
3U0N 
AUON 
3UOQN 
SUON 
3UON 
* i903 


SUON 
3UUON 
3UON 
AUQN 
4UON 
3Uu0N 
3UUN 
4UUN 
SUON 
3UON 
= lod 


SUON 
3SUON 
dUO0N 
3UON 
3UON 
3u0N 
3UON 
SUCN 
3UUN 
3UON 
= Vod4 


3U0N 
au0N 
su0N 
3U0N 
3U0ON 


LS 


an er reer manent 


0°O = AQtoOtazA}ZUSDN_ ATOZ sqtnsey “OT eanhbtd 


© 
We 
& 
© 
PHOOHHHHOHOSHOOHHORAOH OOOO OOD 
(S2s*t- dest" ELlv0" SHSit= 
J19NV 39804 Wwiul A-30¢e03 INVL1NS34Y X-399409 INVLINS3s 
HPHOKOPHOHAHOHHOHHHROGE HOOD DOD 
RL6%s FERS? RYb ER €LT 
*ve°R EES = Ebr’ Ro bul 
deed €eS°S- O¢s"s- 66 
6se°2 ceucR M¢5°bo 64 
17S*t- 9HE° R= 9es*er 6 
*9S°S=- 99n°S~ 067" R- 6 
£66°%l= poK on 2  -d 6S 
EVE R= ees*t foe 6% 
RSo°9- 66%°S 2d¢*s- 6b 7 
SER° a= Ll6cl Lye 6¢ ei 
test cee’s 2loeH- 6t 
Lee's Gong Syren 6 
Ad ¥Z UvOT Iwill IGON 


savul LavHSs 


000000°0 SALTDIGIN3993 HOOK OOOROe 


t 3 e9bL "EL 0€92°2 B09E*= Z2cRB°T vzl= JUON 
Ss 0 1199°8- NO4e°2 BI6E°- 2HeR eT 6tl= 3UON 
0 tu00" Ren2*2 e60%"= SEER TT Bi T= AVON 
se I 6E91 "2 2Hy2°2 980n*= OfEB eT Z,,= 3UON 
t) Ssvo"” REGS2°2 198E°= 2298°1 YLL= SCON 
2 1 2L6L "I~ 09S2°2 B09E*= 2c6u°y Sil= 3uUW 
0 vouo” R952°? 0l62°- R656" 7. L= 3UON 
0 0000"- 0Ls2°2 LUtre"= 6550°? eLl= 3UON 
ee 1 ort" a- 7\y2°2 L601°= Sist*2 ¢ll= 3UON 
) 0 EZ21S° Ot 09L2°2 0000°U Orler? tll= 3U0N 
bo SYlvd dyl3hu 30d04 AWIOVY TNT 4u 1430 JvIave AWILINi® Zl = Vlog 
2iaq elo. "et oEv2°2 8u9E*= 2crB"T ull® 3UON 
s 0 vEBy Be n0%e"2 816¢°= cheesy Sule AUON 
0 9000°- Rene*2 e60%°- Steel vut= 3UON 
€e a 6991 °2 2Be2°2 ROM = 95F8°] dy1= 3UUN 
0 2020" RESC*2 198k °= L296"] yul= SUON 
5 ee 86\8°FI- 09S2°2 RO9E T= Zee T sulle 4UON 
0 0000" 9995¢°2 ol62°- ReS6*y cule JUON 
0 0000" olse%2 (ttZ5> 500°? tul= 3YON 
te I 9SBL"2- 2192°2 lant ’- Siet°? cul= SU0N 
8 0 OvTS*OT NeL2e°2 no000"0 0°12°2 101= 3JUON 
Sore Salva dv 1H3A0 WxHod AWI1IVH TINT 4x 1430 Wwidve AWTLINIA ll = Vos 


wre CS RS . akties ‘RRB RES ; 


*3* 


beto jA) 


i aa :] 


“3°86 


o sgoos: 


sul¥q dV 1HaA0 
S 
0 


t 


0 
6 
) 
9 
0 
0 
dt 
0 

A 


Sdlvd dv 1u3aAV 


- 


onoomoroos 
- 


SuIVd dv 1H3AU 


~ 


omoonomoomnm 
- 


SulVg dV IH3A0 


onoonOMoon 
om 


Sulvd dv 1u3A0 


Z200°O = AqtoTszqUeDDg TOF SITNSOU 


vor2 +17 QLo¢d 2 
0000° oBG2°*2 
0000 "= SHS2°2 
otu6°2> 1292°2 
9leltt NvL2°2 
Wgud WIGVY VWNT da 
16827, Q9e2°2 
6E 1b 6- Si2°2 
€000°> Ty2°2 
Olu6"2 1692°2 
6100°* 99oe"2 
Brbuotm 19s2°2 
oov0'> elsec2 
0000° oL52°2 
Loot 2- +192°2 
LLIS“Ot OyL2°2 
Qodu4d TwIGve WANT du 
0560741 l€92°2 
L9L8°R- Lover2 
T1000" FEt2°2 
Luot2 >B92°2 
8100" Les2°2 
0628" €I~ 9g92°2 
v000°- 1992°2 
00007 £9S2°2 
9€u9 Io Lug92°2 
99U6°6 OvL2°2 
Jodud WIdGVe VWNT da 
l628°el 9272°2 
z2us2"8- 96€2°2 
10v0™ S2r2°2 
Sur} LLy2°2 
6luv "> R2se°2 
v9%S°ETA 9ng2°2 
00v0" 0992°2 
0000" €SS2°2 
4098°= To9et2 
GBS"6 OyL2°2 
30404 WWIdvy WNIT du 
09%S°EL 9l92°2 
60L9"L= 1gEe2°2 
0000° Lty2t2 
7098" tLye*2 
ov00"> 22oe"e 
L2ee*Ft- 19S2°2 
0000° 99S2°2 
0000" ReGe°2 
Bli9°=- R6Sc°2 
8262°4 OvL2°2 
30404 Iwldve WNT du 


1430 


1430 


14330 


1430 


1430 


*TT eaznbtg 


oese - 
GHROf = 
S2l2 
90Tt*= 
o0n0°u 

WIdVe 


72] a 
ecb = 
Sul**= 
S606" 
69RE f= 
Si9e"=- 
Sl6e "=" 
Sites= 
6601°- 
9000° 
VWi1dve 


SIIES= 
S26E"- 
RO0%"= 
BB80%"=- 
o98E°= 
09E° = 
£962" 
eot2°- 
Z60T > 
0000°9 
Widve 


su9E° = 
"l6E"= 
680% >= 
0807" = 
1S8E°= 
%6SE°~ 
2362" 
b602°> 
9R0T*= 
0000°9 
TWidva 


%6SE°- 
So6ee= 
280% = 
olor*- 
GPE o= 
68SE*= 
in62°- 
6R02°~ 
EROT= 
0000°0 
Awildvs 


getned 
HOSOTY 
65°? 
sisi’? 
Ovlc*2 

WWILINT 


zZerB tT 
Zoey] 
Sees eT 
Joey 
LLg8 eT 
2568" 1 
RESET 
6seu°2 
Sisl*2 
Ovi2*2 
AwILInie 


Ze2re"y 
2oreey 
SEE I 
96Eb°1 
LL98°1 
2568 °T 
R6S6"T 
6S9u°2 
Sist*2 
O12°2 
WILINTe 


Ze88"T 
2uoe"] 
Seeh"y 
S6FR TT 
LLOB ] 
eset 
RHS6*T 
6$%0°? 
Sstsic? 
Urle°2 
WILINI+ 


2erkt 
enor? 
Seekel 
FER T 
Ligs*t 
2o68"T 
ReSb TT 
6Sou'? 
sicl*2 
Ove"? 
AWILINIS 


vz00° 


+ 


-awr aye nrnage 
nn 


t 


ALI Ql o4n99923 


4UUN 
3U0N 
JUON 
du0N 


= Mos 


3UON 
3UON 
SUON 
3UUN 
34U0N 
SU0N 
3guN 
3U0N 
3U0ON 
3UON 
= Vos 


3UON 
AUON 
AUON 
SUON 
3UON 
SsU0N 
3U0N 
3U0N 
3UON 
SUON 
= 3108 


3UON 
SUON 
3UON 
3UON 
3UON 
3UON 
AUGN 
3U0N 
S3UUN 
SUu0N 
= 1104 


SUON 
3UO0N 
SUON 
3U0N 
3UUN 
sU0N 
4dUUON 
3U0ON 
3SUUON 
3UON 
= ifo3 


NOJid1eVS aU W4LSAS 


iss 


04 2t 


ree 


LS 


viersapensyiss— 


phn are enna sb 


*7°a 


es e :3 


*9°9 


*°8 


mas el, 


ol 
0 


0 
62 
0 
Salva dv 1N3A0 


- 


n 


oroosornooso 
Nw 


SYuIVd dV 1H3A0 


eoomonrco 


Se 


J 


SYIVd dy TH3A0 


a 


omeoecromoomn 
nw 


SdlVd dv 1d3A0 


Nw 


enrco oomocor 
nN 


SH1Vaq d¥ 1H AAU 


9 
0 
0 
12 
i) 


ee 


Z00°O = AqtTotTazqzUEcO, TOF SzTNSoY 


Qa Steer 
ylbstl- 
1000" 

(@] U8Se"T 
2020° 

I 2915 E17 
0000° 
0000" 

I zZ6u2"2- 
Zsts*ot 

90404 Iwl0Vvy 


Qa TEIs*FL 
SSIT°R- 
1000" 

Qa g802°2 
€020° 

I y6LL°EI- 
Tov0°= 
0000° 

I lu96"2- 
€Zel*tt 

30404 WI10VY 


a SOLL “EL 
6991 B= 
1000° 

a 1096°2 
9500" 

I 8B8E0°+I7 
0000" 
0000°~ 

I oBly°E- 
96LS° Tl 

30404 Ivl0VY 


Q SLE Hl 
1€82°6- 
€000°= 

Q cBle'e 
9500" 

I 96n2°oI- 
0000° 
0000°- 

I 6E99°E- 
Ulett 

3odus Iidve 


(¢] L6H2° 41 
SlE9"6- 
€0vu0'> 

0 0999°° 
9S500° 

I Seve “oIm 
1v00 "> 
0000°> 

I LeSy°e- 
6lis* tt 

Wud IW1dvy 


0 VEE oT 
oLHg*46= 
£000°> 

) SEGH EF 
9500" 


wive?2 
cRE2°2 
91922 
Flo2°2 
1éS2°2 
€SS2°2 
29S2°2 
19S2°2 
o192"2 
OvL2°2 
WNIT de 


f292°2 
e6¢2°2 
F2r2°2 
0842°2 
ors2°2 
“9° 2°e 
glo2°2 
Rlc2°2 
R1I92°2 
OvLe*2 


TNIda 


ve’? 
60%2°2 
2€%2°2 
RBH2°2 
69G2°2 
ls2°e2 
S8S2°2 
RBSG2°% 
4292°2 
OvLz°?e 
WNI4Y 


boye’2 
ely2°2 
Ovy2°2 
b6n2°2 
$952°2 
08s2°2 
T6S2°2 
F6G2°2 
£292°2 
NoL2°2 
JWNIdd 


nSr2°2 
61%2°2 
aone’2 
L6y2°2 
9Gu¢°2 
8s2°2 
T652°2 
7692°2 
9292°2 
Nei2%2 
VWNI dy 


1Sy2°2 
02%2°2 
Grn2°e 
9642"? 
€S42°2 


Se ama” 


“tT eanbty 


26sc °= 
Eooe* = 
ThO# = 
Llov’- 
oSBEe= 
109€°- 
4962°- 
Lot2°- 
S601 "> 
0000*v 
1330 Iidve 


109° = 
l16e"= 
BB0%"= 
980%" 
E9RES = 
219¢"- 
LL62° = 
61 te" = 
eull?- 
0000° 
1430 Wv1dvy 


ZtgEce 
226 "= 
L60%°> 
Z60%°> 
else" 
z29e°> 
LA62°- 
R212°- 
60Tl’> 
0000°0 
1430 WIdves 


Cll 
0000" 


1430 Wy1dve 


B29E°= 
LEese’> 
bole" 
OOtr’= 
6LRE eM 
629E°= 
€662°" 
2E12°= 
110k = 
000070 
1990 WI0ve 


6298 °- 
BE6E = 
bol’ = 
Ab0e° = 
ihe T= 


Zebu 
Zoo8 yy 
Steet 
Goevry 
LL93°1 
25681 
HLS6*] 
6Se0°2 
Sist*2 
Orler? 
WILINITe 


ZZRe"l 
2ote8°T 
Steet 
GhEBT 
Lig8*T 
2s68"l 
ROSE*T 
640°? 
Stsl*2 
Or12°2 
VILINI¢® 


Zcbe"l 
ares] 
SEER *T 
96E8 "I 
Lige Tt 
2cset 
R656" 
6oou"? 
sist*2 
Ovi2"2 
AILINIX 


AITLINIS 


22Pert 
eros t 
Stent 
goeR tT 
Lige't 
2oce°y 
BASb y 
6onu'2 
sisi’? 
Orie’? 
WILINIT= 


arkbet 
dah oe | 
seer 
ater t 
LiSb°] 


vul= 
66 = 
v6 = 
46 = 
yo = 
6 = 
76 = 
t6 
c6 
l6 = 


v6 = 
oe = 
vB = 
ie = 
98 = 
oh = 
oB = 
cn = 
en = 
ig = 


vb = 
ol = 
=! Bie 
it = 
SL = 
ot = 
cL = 
ti = 
ci = 
lz = 


AauON 
sUON 
SUON 
3U0”W 
3uON 
3U0N 
3U0N 
3uoN 
AYON 
3UON 
= W904 


3UON 
3UON 
3U0N 
3UON 
3UUN 
3UUN 
3UON 
3UUN 
3UON 
3UON 
= 1103 


3U0ON 
3UON 
3UON 
3UON 
3UON 
3UON 
SUON 
3UON 
SUON 
3UON 
= W104 


SUUN 
3UON 
S3UUN 
3UUN 
SUON 
4U0N 
3UON 
3U0N 
300N 
AUON 
= od 


3UUN 
SUUN 
SUON 
3U0N 
SUON 
4U0N 
3UUN 
SUUN 
3UON 
3Uu0N 


= liogs 


3U0N 
SON 
30” 
SU0N 
3U0N 


16 


t 
3 
% 


en ae eines 


o 
Z00°O = AQTOTz34USNN0g ZOF SRATNSOY aaa eanbt gs up 
& 
is 3 
ee ee ed 
oozt°e- yoEe cL Llce*d- eesit- 
319Nv 39403 Wwi0l A-306043 LNVLINS3H = X=-39d04 INVL1NSAS 
Cee ee ee ee ee 2a 
168°9 2sn°2 mierda 6Ll 
o9ted 9IE°I- 6letl- 6ul 
LLL°S R16%b= Les*.- 66 
lel °2 OnO*L= 91L"B- 68 
18S*I- £09°8- Lydten ol 
R10°9- R9UL= F8c"6o 69 
0g0°6- lee*e- PLI* OM 6% 
82S°*6- US2*T lsv°6- 6” 
R9T°L= £09 ole 6E 
916°2- £9€°R Liatr= 6¢ Sais 
16%°T rae Oa) ose*6- 61 a 
tl6o°? The’s t29°d- 6 
As x3 uvo7 twill 4Q0N 
sQvUl L3vHSs 
000200° SA) TOTYIN3I09 POORER DOOD 
am! 8zeeel Tl92°2 6HSE°= Zee tT v¢el= 3UON 
s 0 goiecl- TeEee*2 66RE°= 2ceob"l oOLl= 3UON 
0 T000° 2ly2°2 Lloo’- SeEBTT bLl= 3U0N 
SE I RIL9° RIMS 1L09°- 9beerl ZUl= 3U0N 
0 SS00° 12s2°2 bnREen Liseey YLL= SU0ON 
2t 1 96%2°FI- orse°2 ARSE oA 2chel Sil] SUUN 
0 0000° G9S2°2 Lo62°- AEST *L1= 3U0N 
0 0000°- 64s2°2 0602°~ 6S00°2 etl= 3su0N 
€€ I BEBA= 66S2°2 bROT*O sisi’? e1l= 3Y0N 
g 9 9250" ObL2°? 0000°%U Nvicg*? (il= SUON 
"O°8 = SuIV¥g dv 1H3A0 3040s Wave WNT 4x 1430 WIdve WwILIvie el = los 
210 LoH2°ET Nle2°e BHSE TO 2zee"t vil= su0N 
Ss 0 v6le°L-~ oRE2*2 RoRE oH 2reHey by = SUN 
0 9u00°- Ptr2%2 LlLoo'- Steer rule sUON 
ee 0 BERR” AIe7e°2 2lov’- QALH TT dye 3UON 
0 1020" 2g"? Loket= Lissty Yule SUON 
wrt Qree etm #9gz°2 26Se°- aster sul= SUON 
0 0000" 1952°2 €S62°- arheb*y vue au0N 
0 0000° 9S42°2 9602°- bonyt? tul= SUON 
ve Gece" T- Foye"2 RROT"= Sist*2 cul= sUON 
4 0 HLUB TS NoLe°2 0000” Order? tOl= 300% 


bie bel) Sulvd dv 143A0 39404 TWlOVeH TNT 44 1330 TW1dvo VWTLINIt 11 = Wos 


components are presented. The resultant forces acting on the 
Shaft are also specified, i.e., EF Ss cP ye total resulting force 
and attitude "angle". 


2.7 Input Data Presentation (Figures 14 and 15) 


Card I - Title Card Format (10A8) 


Provided for a descriptive title of the problem being 


analyzed. 
Card II - Foil Geometry Card Format (8F10.5) 
Word 1 - Number of foils 


Word 2 - Number of nodes per foil (6, 8 or 10) 


Word 3 - Radius of bearing (inches) 


Word 4 — Axial length of bearing (inches) 


Word 5 - Thickness of elastic foil base metal (inches) 


Word 6 - Thickness of nonstructural coating (inches) 


SRT 


Word 7 - Arc length of foil (inches) 


= 


Word 8 - Free foil radius (inches) 
Card III - Plotting Card for Future Program Development 


"Leave Blank" 
Card IV - Foil Modulus of Elasticity Card Format (F10.0) 


Word 1 = Modulus of elasticity (psi x 10°°, 
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Card V - Shaft Radius Card Format (F10.0) 
Word 1 ~ Shaft radius (inches) 
Card VI - Eccentricity Control Card Format (3F10.0) 


Word 1 - Final eccentricity (inches) 


Word 2 - Increment of eccentricity (inches) 
Word 3 - Initial eccentricity (inches) 
Card VII - Attitude Control Card Format (F10.0) 


Word 1 - Angle (¥) of Figure 5 (Radians) 


2.8 Comparison of Analytical and Experimental Data 


Numerous test cases have been evaluated to demonstrate the 
programs capability in recognizing proper release, kinematic con- 
straints, and foil overlap in duplicating experimental static 
stiffness data. A thorough examination of the elasticity program 


and correlation with experimental data has demonstrated the ade- 


quacy of this portion of the elasto-hydrodynamics journal bearing 
program. Some of these results are presented below: 


Se PT aaa 


B fe) Figure 16 illustrates a typical foil geometry and 
Ps the coordinate system utilized in this study. Note 
that the positive direction of eccentricity (e) has 
been established in the direction of the positive Y 
axis. Accordingly, a negative value of (Fy) on the 
journal opposes a positive value of eccentricity. 
Further, a positive value of (FY) encourages forward 
precession (whirl), whereas a negative value of 
(FL) opposes forward whirl. The angle (¥) is a 
geometric parameter and is a user input option to 
control the displacement vector (eccentricity) 


direction with respect to foil geometry. 
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Figure 16. Geometry and Coordinate System | 
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Figure 17 shows a comparison of experimental data 
and analytical data for a 0.0085 inch thick 16- 


segment foil coated with 0.002 inch of Teflon-S. 


As shown, excellent correlation was achieved. 


Figure 18 is a comparison of analytic and 
experimental data for a 0.010 thick uncoated 8- 
segment foil. Two curves are present; one for 

the displacement vector (eccentricity) directed 
toward the trailing edge (T.E.) of a foil, the 
other with the eccentricity directed midway between 
two trailing edges of adjacent foils. It is inter- 
esting that both experimental and analytic data 
demonstrate the "stiffer" curve is produced when 
loading is directed midway between the trailing 


edges of two adjacent foils. 


Figure 19 presents the data obtained for a 0.012 
inch thick 12-segment foil coated with 0.001 inch 
of Teflon-S. Two curves for the experimental and 
analytic study are presented similar to those of 
Figure 18 for the 8-segment foil. It is noteworthy 
that unlike Figure 18, the experimental curves for 
loading over the trailing edge and loading between 
adjacent foils cross at approximately 0.005 inch 
eccentricity. This same characteristic is dupli- 
cated at an eccentricity of 0.006 inch for the 


analytic data. 


Figure 20 presents the loading obtained analytically, 


which represents the resultant orthogonal load FS. 
(see Figure 5), for the 12-segment foil studied in 
Figure 19. The magnitude of this load is small in 


comparison to the load-resisting eccentricity (FY) 


of Figure 19. It is interesting to note the apparent 


symmetry when the loading direction is altered. 
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Figure 21 is a 12-segment design analytic study of 
the effect of foil length on static stiffness. The 
solid line represents a minimum foil length (0.010 
inch less than nominal) or a foil overlap of 47.2 
percent. The "stiffer" dashed curve illustrates 

the load versus deflection data for a maximum foil 
length (0.010 inch greater than nominal) or an over- 
lap of 47.7 percent (the nominal foil overlap is 


47.5 percent). 
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SECTION III 


JOURNAL BEARING ANALYTIC PROGRAM 


Sei General 


The foil journal bearing computer elasto-hydrodynamic analy- 
Sis parallels the thrust bearing computer analysis program. A 
separate hydrodynamics subroutine and elasticity subroutine ex- 
change the appropriate information on loads and foil deformations 
until convergence is achieved. Convergence is attained when the 
hydrodynamic film thickness used to obtain the pressure distribu- 
tion and resulting foil loading is the same film thickness (with- 
in a reasonable error) predicted by the elasticity program when 
the hydrodynamic loads are imposed on the foils. The Foil Journal 
Bearing Program flow chart is provided in Figure 22. One funda- 
mental but significant feature of journal bearing geometry, is 
that journal foil segments are not independent elastic members 
but are in contact with each other. Thus, each foil encounters 
loading due to the hydrodynamics of the bearing as well as trac- 
tion loads and kinematic constraints from foil interaction. [In 
simple terms, foils may: (1) be in contact with each other, (2) 
be in contact with each other and constrained by the bearing hous- 
ing, (3) be in contact with boundaries and not in contact with 
other foils, or (4) be in contact neither with other foils nor 
bwandaries. These conditions must be correctly assessed by the 
elasticity program, thereby requiring a difficult nonlinear iter- 
ative elasticity subroutine to be programmed. This has been 


successfully achieved. 


3.2 Discussion of Geometry Subroutine 


As all foil displacements must be referenced from pre- 
scribed datum, it is first necessary to determine this reference 
configuration. For the foil thrust bearing this datum is the 


free foil geometry. For the journal bearing, the configuration 
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“+3.3 Elasticity Subroutine 


of the foils may differ in the assembled state from the free form 
state due to foil interaction loading. Generally, the fewer num- 
ber of segments, the more likely that the assembly configuration 
deviates from the free foil state. The assembly subroutine has 
been programmed to determine the assembled configuration of the 
foil journal bearing prior to shaft insertion. This geometry 
subroutine has been programmed to search for an acceptable con- 
stant curvature position at assembly. Given the free foil radius 
(RFF), the length of the foil (FL), the base material thickness | 
(te), the non-structural coating thickness (tis and the bearing 
radius (Rp) the assembly subroutine seeks a foil radius Ro (Fig- 
ure 23) where R' = R + te + te The assembly routine reports Rot 
as well as Ry (the radius to each node from the geometric center 
of the bearing) from which all foil deformations are measured. 
Figure 24 presents some typical Calcomp computer plots generated 


by the assembly routine for 6-, 12- and 16-segment foils. 


Figure 25 presents the details of the assembly/elasticity 
solution. As noted, the elasticity solution commences by obtain- 
ing the symmetric "preload" position of the foils. This is the 
configuration when the shaft is inserted in a zero "g" position. 
The program proceeds from there to obtain the foil configuration 
when the hydrodynamic loads are applied. Program checks provide 


for: 


(a) Constraint or release of nodes on adjacent foils 


and/or contact with boundaries 


(b) Constraint or release of nodes in contact with 
the shaft 


(c) Constraint or release of nodes on the bearing 


housing 
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Nodal constraint is provided when displacements indicate a 
violation of kinematic boundaries. Release of nodes occurs when 
traction loads produced by constraint are no longer physically 


realizable. 


The foil elasticity is analytically represented as a series 
of short curved elements with a coordinate system and sign con- 


vention (Figure 26). 


By application of Castigliano's second theorem, the element 
and system stiffness matrices may be determined. The stiffness 
equation solution is provided by a straight-forward displacement 
matrix approach. The attachment of the foil to the bearing hous-~ 
ing is assumed to be a "pinned joint". The difficulty in the 
solution is encountered when the foils are loaded. The foils 
begin to contact each other and encounter kinematic constraints 
of the journal and bearing housing. It is apparent that this 
provides a highly nonlinear stiffness characteristic of foil 
journal bearings. To deal with these problems the elasticity 
program generates a constraint matrix. This constraint matrix is 
used in a transformation of the governing equations to deflate 
the set of equations. By proper partitioning of this reduced set 
and by back transformation, the solution to the constrained equa- 
tion is determined. This permits the nonlinear stiffness char- 


acteristics to be simulated. 


J3e4 Hydrodynamic Analyses 


The journal foil bearing subroutine termed the “Hydrodynamics 
Program" (Figure 22) is a finite difference solution to the 


Reynolds equaticn, based upon the assumptions that: 


(a) The film clearance between the shaft and foil is 


small compared to the radius and length of the bearing. 


x a2500 


a =e 


mess £5 


<« 


i 


NODE N+1l 


NODE N 


= Radius of curvature 

= Tangential load 

= Bending moment 

= Shear load 

= Angle which subtends the element 
= Tangential displacement 

= Flexure displacement 

= Rotation 


cos &£ Ww OO & VU DW 
I 


= Thickness 


Figure 26.. Foil Geometry Coordinate System. 
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(b) 


(c) 


(a) 


(e) 


(f) 


(g) 


(h) 


the hydrodynamics solution. 


and radial velocity components U and V, 


The flow may be laminar or turbulent depending upon 


the local Reynolds number. 


Wet surfaces - no slip at bearing surfaces is 


permitted. 


The working fluid may be compressible or 


incompressible. 


The velocity gradients across the fluid film are 


larger than other velocity gradients. 


Axial and circumferential flows are permitted. 


Program solution for finite length bearings. 


Body forces, viscosity and density may vary with 


position in the axial and tangential directions. 


Nonsynchronous and squeeze film loading permitted. 


Figure 27 illustrates the foil journal bearing geometry for 


respectively. 


pressed as: 


' 3 3 
sf» n° foe . = a {fe a pee fel. 
ale (3 )|+s iE 7 Me 2) 
vA 
G6 (PG. -,) sh + Ph s (UL, + UL) + Rh (U, + 0, )8P 
(PU, 1’ %& a ‘i ed) 
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Figure 28 illustrates the tangential 
With this 


geometry defined, the compressible Reynolds equation can be ex- 


ee 


‘ x',y' FIXED FRAME AXES 
x, yY WHIRL FRAME AXES 


Figure 27. Foil Journal Bearing Geometry for 
Hydrodynamic Solution. | 
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Figure 28. Foil Journal Bearing Tangential and 
Radial Velocity Components. 


ee. 


where: 
P = hydrodynamic pressure 
0 = Density = P/RT (Isothermal conditions) 
X, Z = body forces 
Use V5 = velocity components at shaft surface 
Uy, vy = velocity components at foil surface 


h = film thickness 


yu = viscosity 


and the turbulence constants k, and ko are given by: 


z, e 2 0.725 
k. = 12 + 0.53°(K R,) 


5 


Ul 


ee 0.6 
k L2 + 0.296(K R 
‘ (K°R,) 


L ~ where Ra is the local Reynolds Number and K is the Prandtl's mix- 
Pe ing length. 
: 3.5 Comparison of Analytical and Experimental Data for 


Elasticity Solution 


Numerous test cases have been evaluated to demonstrate the 
program's capability in recognizing proper release, kinematic 
constraints, and foil overlap in duplicating experimental static 
stiffness data. A thorough examination of the elasticity program 
and correlation with experimental data has demonstrated the ade- 
quacy of this portion of the elasto-hydrodynamics journal bearing 


program. Some of these results are presented below: 
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fe) Figure 29 illustrates a typical foil geometry and 


the coordinate system utilized in this study. Note 
that the positive direction of eccentricity (e) has 
been established in the direction of the positive Y 
axis. Accordingly, a negative value of (FY) on the 
journal opposes a positive value of eccentricity (e). 
Further, a positive value of ae encourages forward 
precession (whirl), whereas a negative value of (PL) ] 
Opposes forward whirl. The angle (¥) is a geometric 

parameter and is a user input option to control the 

displacement vector (eccentricity) direction with 


respect to foil geometry. 


fe) Figure 30 shows a comparison of experimental 
data and analytic data for a 0.0085-in. thick 
16-segment foil coated with 0.001 in. of Teflon-s. 


As shown, excellent correlation was achieved. 


° Figure 31 is a comparison of analytic and ex- 
perimental data for a 0.0i0-in. thick uncoated 
8-segment foil. Two curves are presented; one 
for the displacement vector (eccentricity) dir- 
ected toward the trailing edge (T.E.) of a foil, 
the other with the eccentricity directed midway 


between two trailing edges of adjacent foils. It 


erg: PFs 


is interesting that both experimental and analytic 
data demonstrate the "stiffer" curve is produced 
when loading is directed midway between the trail- 


ing edges of two adjacent foils. 


fe) Figure 32 presents the data obtained for a 0.012-in. 
thick 12-segment foil coated with 0.001 in. of 
Teflon-S. Two curves for the experimental and an- 
alytic study are presented similar to those of 
Figure 31 for the 8-segment foil. It is noteworthy 


that, unlike Figure 3l, the experimental curves for 
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Figure 29. Geometry and Coordinate System for Elasticity 
Program 
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am!) 


loading over the trailing edge and loading between 
adjacent foils cross at approximately 0.005-in, 


eccentricity. This same characteristic is duplicated 


at an eccentricity of 0.006 in. for the analytic data. 


fe) Figure 33 presents the loading obtained analytically, 
which represents the resultant orthogonal load eS 
(see Figure 29), for the 12-segment foil studied in 
Figure 32. The magnitude of this load is small in 
comparison to the load-resisting eccentricity (FY) 
of Figure 32. It is interesting to note the apparent 


symmetry when the loading direction is altered. 


fe) Figure 34 is a 12-segment design analytic study of 
| the effect of foil length on static stiffness. The 
solid line represents a minimum foil length (0.010 
| in. less than nominal) or a foil overlap of 47.2 per- 
cent. The “stiffer” dashed curve illustrates the load 
versus deflection data for a maximum foil length (0.010 
F in. greater than nominal) or an overlap of 47.7 percent 


(the nominal foil overlap is 47.5 percent). 


: 3.6 Foil Bearing Analytical Program Preliminary Results 
Preliminary results of the nonlinear elasto-hydrodynamics 
program provided solutions that converged to an acceptable 
elasticity configuration, as well as an acceptable hydrodynamics 


solution. 


Figure 35 represents the symmetric pressure profile for the 
zero eccentricity case of an 8~segment foil bearing. The com- 
4 puter solution verified the no~load configuration and the symmetry 


of pressure distribution. 
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Figure 36 is the same 8-segment configuration with an eccen- 


tricity of 0.004 in. The nonsymmetric pressure profile illus- 


trates the pressure loading to oppose eccentricity. 


3.7 Journal Bearing Computer Program 


This paragraph presents the input data code foliowed by a 
sample configuration that has been analyzed by the CDC-6400 Pro- 


gram. A Fortran program listing is provided in Appendix B. 
3.7.1 Input Data Presentation 
Card 1 - Title Card ~ Elasticity Format (8Al10) 


Provided for a descriptive title of the problem being 


analyzed. 

Card 2 - Foil Geometry Card Format (8F10.5) 
Word 1 = Number of foils 
Word 2 - Number of nodes per foil (6, 8 or 10) 
Word 3 - Radius of bearing (inches) 
Word 4 - Axial length of bearing (inches) 
Word 5 - Thickness of elastic foil base metal (inches) 
Word 6 - Thickness of nonstructural coating (inches) 
Word 7 - Total arc length of foil (inches) 


Word 8 - Free foil radius (inches) 


FILM 


Figure 36. 


PRESSURE PROFILE 


2.250 
0.010 IN. 
0.001 IN. 
2.282 IN. 
0.004 IN. 
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MINIMUM FILM 
0.00007 IN. 
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3.36 LB-IN. 


Eight-Segment Foil Center Plane Pressure Distribution 
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Card 


Card 


Card 


Card 


Card 


Card 


Card 


3 - Plotting Card for Future Program Development 
"Leave Blank" 

4 - Modulus of Elasticity Card Format (F10.0) 
Word 1 - Foil modulus of elasticity (psi x 10°°) 

5 - Shaft Radius Card Format (F10.0) 
Word 1 - Shaft radius (inches) 

6 - Eccentricity Control Card Format (3F10.0) 
Word 1 - Final eccentricity (inches) 

Word 2 - Increment of eccentricity (inches) 

Word 3 - Initial eccentricity (inches) 

7 - Attitude Control Card Format (F10.0) 
Word 1 - Angle (¥) of Figure 29 (Radians) 

8 - Title Card - Hydrodynamics Program Format (8A10) 

9 - Fluid Definition Card Format (6E12.3) 
Word 1 - Fluid Density (1bs/in.°?) 


Word 2 Fluid Viscosity (1bs-sec/in.*) 


Word 3 - Leave blank for foil bearing solution 


Word 4 - Relaxation Factor for Hydrodynamics - 
Elasticity Iterations (Recommend 0.25) 
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Card 


Card 


Word 


Word 


Word 


Word 


Word 


Word 


Word 


NOU 


Word 


Word 


Word 


Word 


Word 


Word 


5 - Permissible Convergence Error for Solution 


(Recommend 0.1 to 0.25) 


Boundary Conditions Card Format (5E12.3, F6.0) 


1 - Boundary Pressure Row 1 (psi) 
(Rows made by longitudinal lines parallel 
to centerline) 

2 - Boundary Pressure Last Row (psi) 


3 - Vapor Pressure of Fluid (psi) 


4 - Error Permissible for Pressure Iterations 


of Hydrodynamic Solution (Recommend 1.0) 


5 - Factor to Over-or-Under-Relax Pressure 


Iterations (Recommend 1.0) 


6 - Number of Iterations for Extrapolation 


(Recommend 5.0) 


Bearing Kinematics Card Format (6E12.3) 

1 - Spin Speed of Shaft (rpm) 

2 - Whirl Fraction (Multiplied by Spin Speed to 
Obtain Whirl speed ~- a value of 1.0 indicates 
synchronous whirl) 

3 - Leave Blank 

4 - Time Derivative of Eccentricity (in/sec) 


5 - Value of a (Figure 27) (Degrees) 


6 - Time Derivative of a (rpm) 
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i Card 12 - 


Word 
Word 
4 Word 
2 Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 
Word 


Card 13 - 


M Word 
Word 
Py Word 
” Word 


no 


Word 


To illustrate the use of this computer code, 


input and results for the geometry depicted 


Hydrodynamic Control Card 


wo ON DO FW DY KF 


el ee ee ee ee 
on OO WwW FP WD FF OO 


Geometry Card 


& W DY Ff 


Format (1814) 


Number of Axial Rows 


Number of Circumferential Columns 


Leave Blank 


= 500 


Leave Blank 


=1 


Il 
WDM FP HN FN DY NY W F&F ND 


Row Number 


These values must be specified 
as indicated for proper oper- 
ation of the Hydrodynamic 
Solution for Foil Bearing 


analysis 


Format (215, E12.3, 15, I2) 


(I) 


Column Number (J) 


Specified Pressure at (I, J) 

Code 

0 = Interior Point 

1 = Boundary Point Number Iterated 
Next Card Code 

1 = Continue reading 

2 = Last Card 


a ee: Sample Problem 


presented on the following pages. 
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in Figure 37, 


the computer 


are 


SLT 


No. of Elasticity Nodes = 80 
Hydrodynamic Grid = 35 x 5 

No. of Foils = 8 

No. of Elasticity Nodes/Foil = 10 
Bearing Radius = 2.282 in. 

Axial Length = 6.000 in. 

Foil Thickness = 0.010 in. 

Coating Thickness = 0.001 in. 

Arc Length of Foil = 3.40 in. 

Free Form Radius = 3.00 in. 

P = 15 psi (Ambient) 

Eccentricity = 0.004 in. 

Spin Speed = 33,000 rpm 

Viscosity = 0.42 x io ° lbs/sec/in? 
Density of Fluid = 2.29 x ios lbs/in? 
Modulus of Elasticity = 30 x ae” psi 

Y = 0.0 Radians 


Figure 37. Sample Problem Geometry 
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SECTION IV 


THRUST BEARING ANALYTIC PROGRAM 


4.1 General 


A coupled elasticity and hydrodynamics computer simulation 
for the self acting compliant foil thrust bearing has been com- 
pleted. This analytic program has been experimentally verified 
and has proved to be a valuable aid in permitting increased load 
capacity to be achieved. This program represents a major first- 


step in foil thrust bearing analytic capabilities. 


This section discusses the analytic basis with examples of 
input and program results. In addition, a detailes discussion | 


of analysis and experimental results is provided. 
4.2 Thrust Bearing Analytical Program 


The Thrust Foil Bearing Analytic Solution consists of a 
finite difference hydrodynamic analysis and a finite element 
elasticity routine, which exchange information to obtain a con- 
verged solution. As shown in the program logic diagram, Figure 
38, the program commences by requiring a definition of the bear- 
ing geometry and initial film thickness. The starting film 
thickness is a user input requirement. The closer the initial 
film is to the converged solution, the more rapid the convergence 
of the elastohydrodynamics solution. Once the initial film and 
input parameters are defined, the program proceeds to the hydro- 
dynamic analysis to determine a fluid pressure distribution and 
the resulting loads to be applied to the finite element elasti- 
city routine. With these loads, the elasticity solution provides 
the foil deflections and consequently, the elasticity film dis- 
tribution (h,)- With (h) and the hydrodynamic film (hy) de- 


fined, a relaxation of the form: 


101 


*wearBbetq otboy weahoag soOtweuAporzpAH-o7selTqA “ge oanbtgy 


SSANAYOTHL 
Wild IWNId 


ee se WOaHO 
r L__Aowsnoov 


MOWHO 
AONVI'TdWOO 


NOTLOaTAdd 
TI0O4 


X1IaLVW 
SSANIALLS 
' TIO 

wusoud =F —-—---——-— 


ALIOILSV1a qiala AYASsaud 
avOT DNIWVad 


NOT LVOTAIGOW 
NOILIGNOO 
AUVGNNO 


cata! 


WWdooud 
SOINWNAGOUCAH 


. 
| 
| 
| 
{ 
I 

¥v 


SSANMOTHL 
WIId 
TIWILINI 


102 


(1) (i-1) 


3 ay. 
: Ne = ah, + (1 a) hy 


Where: 


(i) is the iteration index 


a@ is the relaxation parameter 


permits redefinition of a new hydrodynamic film and the cycle 

continues. When the film distribution used in the hydrodynamics 
(Dy 
@. ) then convergence is achieved as: 


anh) is the same as that predicted by the elasticity routine 


(i) a aT) a GbR Lt) 
he ee = hy 


4.3 Elasticity Solution 


The elasticity solution is a three-dimensional triangular- 
plate finite element representation of the foil with provisions 
for assigning spring stiffness at prescribed nodal locations. 
This linear elasticity solution contains provisions for 2valua- 
ting thermal effects, including that due to thermal gradients 
across the plate thickness. The elasticity routine accepts the 
pressure distribution, as determined by the hydrodynamic finite 
difference solution. The pressure loading is then integrated for 


nodal loads, which provide foil deformations for the linear elas- 


UP PLE TL 


ticity solution. Special boundary conditions such as symmetry 
are directly provided by proper conditioning of the finite ele- 
ment elasticity matrix. A schematic of a typical thrust bearing 
model, illustrating the finite difference and finite element 


grids is shown in Figure 39. 


The finite element grid consists of triangular shell ele- 
ments. As shown in Figure 40, three cartesian components 
of translation and rotations at each of the three mid-thickness 


vertex are provided. 
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Figure 39. Thrust Foil Bearing Mathematical Model. 
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Figure 40, Finite Element Grid. 

The local coordinate system shown in Figure 40 is utilized 
in the computations but foil elasticity deflections are reported 
by the program in a global coordinate system. 

4.4 Hydrodynamics Solution 

The thrust foil bearing subroutine termed the hydrodynamics 

program of Figure 38 is a finite difference solution to the 


Reynolds equation. Based upon the assumptions that: 


(a) The axial clearance between bearing surfaces is 


much smaller than other bearing dimensions. 


(b) The flow may be laminar or turbulent. 


(c) Both bearing surfaces are wet. 


(d) The velocity gradients across the fluid film are 
4 much larger than other velocity gradients. 


| (e) Body forces, pressures, viscosity and density are 


invariant across the lubricant film. 
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(£) The working fluid may be compressible or 


incompressible. 
(g) o = P/RT (Isothermal conditions). 


The Reynold's equation can be expressed in a cartesian frame 


as: 


Where subscripts (0) and (1) refer to velocity components at the 


surfaces indicated on Figure 41 and: 


P = the pressure 
op = the density 
X, GZ = body forces 
V, = velocity components 
h = film thickness 


yw = viscosity 


and the turbulence constants ky and ko are given by: 


k= 22 + O53 oe 


kK, = 12 + 0.296 waa 
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Figure 41. Thrust Bearing Geometry. 
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Where Ra is the local Reynolds Number and Kappa (K) is the 
Prandtl's mixing length, which has been experimentally deter-~ 
mined to be dependent upon the eccentricity ratio for plain 
journal bearings as shown in Figure 42 with analytic correla~- 


tion per Figure 43. 


Transformed to cylindrical coordinates, and accounting for the 


velocity terms, the Reynold's equation becomes: 


3 poe - ap (Pp ni (ap _ = 
— | B foes Bl + oS le R 
rag ky. u rae ox ku ax 


PU hU_ oP 
6 (=2 oe + 2, 


r P) 2 hoe 
where: 
v= a - rf cos(a-90) - r?g SIN(6-90) 
UL = xrQ 


Q = rotation speed 


@ is the angular location of minimum film 


clearance 
e implies derivative with respect to time 


The Reynolds equation is written and solved in finite 


difference form to provide an expression for pressure at the 


point (1,7). 
L+1,j 
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Pressure distribution within the bearing is obtained by an 
iteration technique. The program includes the option of consid- 
ering fluid inertia. Mean values of the fluid inertia terms ae 
and ope are computed utilizing pressure distributions at inter- 
mediate stages of iteration and incorporated in subsequent itera- 


tions until the designated accuracy requirement is satisfied. 
4.5 Computational Results 


The baseline foil thrust bearing, as specified below, was 
first analyzed using this newly developed elasto-hydrodynamics 
computer program. Figures 44 and 45 show typical computer re- 
sults of film thickness and film pressure distribution in iso- 
pleth form for one thrust bearing sector at a bearing load of 


100 pounds (5 psi). 


Outside Diameter (OD) 6.15 In. 
Inside Diameter (ID) 3.65 En. 
Thickness of Foil (T)) 0.006 In. 
Thickness of Foil Pad Plate (T,) O00) In. 
Thickness of Stiffener (T3) 0.012 In. 
Thickness of Spring Pad (Ty) On Olu2 In. 
Number of Sectors (N) 3 


Figure 46 shows film thickness results and runner positions 
at the mean radius (R = Rw for various bearing loading condi- 
tions. The effects of foil bearing compliance are also shown on 
Figure 46. Since the film shape becomes extremely steep at 
higher loading conditions (Figure 46), higher friction losses are 
expected. The bearing friction loss versus bearing load, as pre-= 
dicted by the analytical program is shown in Figure 47. Gas film 
pressure profiles along the mean radius line (R = R.) for various 
bearing loading conditions are shown in Figure 48. The minimum 
film thickness, usually located near the bearing inside diameter 


(see Figure 44), is shown against the bearing load in Figure 49. 
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Figure 44. Baseline Bearing Film Thickness Isopleth 
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Figure 45. Baseline Bearing Film Pressure Isopleth (psi) 
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Figure 46. Film Thickness at Ru for Baseline Thrust Bearing 
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Figure 49- Baseline Bearing Load as Function of Minimum 
Film Thickness 
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Bearing load capacity is predicted by using a minimum film 


thickness of 0.0001 inch. The baseline bearing was predicted to 


fail at a bearing load of 280 pounds by the analytical study. 


The baseline thrust bearing was tested in the test rig. The 
bearing carried a load of 304 pounds. The rubbing pattern shown 
in that test clearly indicated that rubbing was initiated at the 
regions of minimum film thickness as predicted by the analytical 


results. 
4.6 Thrust Bearing Computer Program 


This paragraph presents the input data code followed by a 
Sample configuration that has been analyzed by the CDC-6400 Pro- 
gram (Figure 50). A Fortran program listing is presented in 


Appendix C. 
4.6.1 Input Data Presentation 


eCard iL Format (10A8) 
Word 1 - Comment card for hydrodynamics subroutine 
Card 2 Format (6E 12.3) 


- Inside radius of bearing (inches) 
- Maximum radius of bearing (inches) 
Fluid density (1bs£/in°) 

- Fluid viscosity (1bsf-sec/in”) 


- Vapor pressure (psi) 


= 
ie) 
n 
Q 
Dn oO Be WN 
1 


- Ambient pressure (psi) 


Card 3 Format (6E 12.3) 


Word 1 - Turbulence constant (A value of 0.2 is recom- 
mended) 


Word 2 - Critical turbulent Reynolds number 
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Card 


Card 


Ww NY F 


ND WU Bf 


16 


Allowable extrapolation pressure error 
(Recommend 0.001) 

Pad span angle (degrees) 

Hydrodynamic relaxation factor 
(Recommend 0.5) 


Minimum number of iterations (Default is 3) 
Format (6E12.3, 6X, 12) 


Rotational speed (rpm) 
Pivot coordinate radius (in) 


Pivot coordinate angle (degrees from leading 
edge to center of pressure) 


Pad radial slope (Geometric, use zero) 
Pad tangential slope 
Pivot clearance (in) 


-Not used- 
Format (1814) 


No. of rows (m) 
No. of columns (n) 
Limit pressure iterations (Recommend 500) 


Print number of iterations for pressure 
(1 = Print, 2 = No print) 


2 

Zz 

1 

2 These values must be specified 

2 as indicated for proper operation 
2 of the hydrodydynamic solution 
2 for foil bearing analysis 
iL 


Read flow (Incompressible = l, 
Compressible = 2) 


Fluid inertia (l = yes, 2 = no) 
Limit inertia iterations = 5 
2 


119 


62500 | 
| 


Card 


Card 


Card 


Card 


Card 


Card 


Card 


Word 17 = 2 
Word 18 = 2 
5(0)* Format (48X, 2F 12.3) 


Optional card for compressibie flow 
Word 1 - Gas constant ft/°R 


Word 2 = Gas temperature °R 


6 Format (10 F7.2) 
Words 1 through n ~ Read a5 through a 
(n defined as Word 2 of Card 5) 
7 Format (205, El2cs, £2, 26) 
(Number of cards as required to describe each nodal 
correspondance between elasticity and hydrodynamics 
nodes. See Figure 50 for a typical grid) 
Word 1 - Hydrodynamic Node n 
Word 2 - Hydrodynamic Node m 
Word 3 - Film clearance of node 
Word 4 - Next card 1 = more cards to be read 


2 = last card 
Word 5 = Elasticity node number 


8 Format (E12.3) 
Word 1 - Maximum permissible film (in) 


9 
EOR - End of record (7/8/9) for hydrodynamics program 


10 Format (10A8) 
"Stiffness" 

bi Format (10A8) 
"Ena" 


*This card is used only if Word 13, Card 5 equals 2. 
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Cards 12 and 13 Format (10A8) 
Comment cards for elasticity selution 

Card 14 Format (10F 7.0) 
Word 1 - Plate modulus of elasticity x 10° psi 
Word 2 - Number of specified springs per foil 
Word 3 = Poisson's ratio of plate 
Word 4 - Density (default = 0.16) lbs/in° 
Word 5 = RPM in X direction 
Word 6 - Skew of X-Y plane to X axis of rotation 
Word 7 - RPM in Z direction 
Word 8 - Coefficient of thermal expansion (in/in °F) 


of plate 
Word 9 —- Relaxation factor for elasto-hydrodynamic 
iteration (recommend 0.12) 
Word 10 - Number of special boundary conditions - 
Displacement nodal pairs (see Card 18) 


Card 15 Format (F7.0, 9F9.0, 1X, D1, 
Tx) 


(One for each elasticity node) 


Word 1 - Elasticity node number 


Word 2 = Ux 

Word 3 = Uy 
Word 4 = Uz Boundary condition code 

Word 5 = 8x 

Word 6 = @y 

Word 7 =- @z 

Word 8 - x 

Word 9 = y Cartesian coordinates of node 


Word 10 - z 
Word 11 - Last card 0 = more cards 15 to follow 
l== last card 
Word 12 - = 1 for skew symmetry independent node 
(triggers dependency of node (Word 1) for 
Card 18) | 
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Card 


Card 


Card 


Card 


16 Format (F7.0, 9F9.0, 1X, I1) 


(One for each triangular elasticity element) 


Word 1 - Element Number 

Word 2 - Node l 

Word 3 - Node 2? £=Three nodes watch for vortices of 
Word & = Node 3 triangular plate 

Word 5 - Element thickness (in.) 

Word 6 - Modulus of elasticity of plate 


(Default word 1. card 14) 
Word 7 - Poisson's ratio of plate (default Word 3 - 
Card 14) 
Word 8 - Temperature of top surface of plate (°F) 
Word 9 - Temperature of bottom surface (°F) 
Word 10 - Coefficient of thermal expansion of plate 
(default Word 8 - Card 14) 
Word 11 - Last card indicator 0 = another card to follow 


1 = last card 


uly Format (7F7.0) 
(Number of cards per Word 2 - Card 14) 

Word 1 - Node with specified spring 

Word 2 - Kx 


Word 3 - Ky Spring stiffness (lbsf/in) 
Word 4 - Kz | 
Word 5 = K@x 
Es Rotational spring stiffness 
Word 9" BOY | Gimelb/tadban) 
Word 7 - K®éz 
18 Format (2F7.0) 


(Number of cards per word 10 - Card 14) 
Word 1 - Node 1| 
Word 2 = Node 2 


Displacement pair nodes 


19 Format (F7.0) 
Word 1 - Number of Cards 20 
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Card 20 Format (6F7.0) 


Word 1 - Ux 
Word 2 - Uy 
Word 3 - Uz Displacement preload on springs 
Word 4 - @x 
Word 5 - 6y 
Word 6 - 902 
Card 21 


EOR - End of Record (7/8/9) for elasticity program 


(One for each elasticity node) 
Word 1 -— Node number 


Word 2 — Initial hydrodynamic film (in) 


4.6.2 Sample Problem 


To illustrate the use of this computer code, the computer 
Pa input and results for the configuration depicted in Figure 50, 
* are presented on the following pages. 
é] 
Ba 
| 
- i 
123 


e25900 


Card 22 Format (1X, 15, F14.8) 


(49) (48) 


ELASTICITY 
ELEMENT NO. 


Pew 


HYDRODYNAMIC NODES n x m 
ELASTICITY NODES ( ) 
SPRING LOCATIONS = 

27.692° 


Figure 50. Sample Thrust Bearing Configuration 
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